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Using a combination of analytic arguments and numerical simulations, we determine lower and 
upper bounds for the energy barriers to the motion of a defect line in a random potential at low 
temperatures. We study the cases of magnetic flux lines in high-T c superconductors in 2 and 3 
dimensions, and of domain walls in 2 dimensional random-field Ising models. The results show that, 
under fairly general conditions, energy barriers have the same scaling as the fluctuations in free 
energy, except for possible logarithmic factors. This holds not only for barriers between optimal 
configurations of the line, but also for barriers separating any metastable configuration from a 
configuration of minimal energy. Similar arguments may be applicable to other elastic media with 
impurities, such as bunches of flux lines. 
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I. INTRODUCTION 

Directed paths in random media (DPRM) Q] are sim- 
ple realizations of glassy systems ^|||. Some examples 
are pinned flux lines (FL) in high-T c superconductors, 
and domain walls in random field and random bond Ising 
models. In thermal equilibrium, a magnetic FL is pinned 
by defects (oxygen impurities, grain boundaries, etc.) in 
the superconductor which lower its energy H. The re- 
sulting elastic distortions are limited by the line tension 
which opposes the bending of the line. This comptetition 
leads to a free energy landscape for the FL which is rather 
complicated and has many local minima, i.e. metastable 
states When an electric current flows through the 
system, the FL feels a Lorentz force perpendicular to 
its orientation and to the current direction. As long as 
the current is not strong enough to overcome the pin- 
ning forces, the line moves by thermally activated jumps 
of line segments between metastable configurations ||-|) . 
The length of these line segments is estimated by the con- 
dition that the free energy barrier for a jump should be 
of the same order as the gain in free energy due to that 
jump. These dynamics are believed to be the reason for 
the nonlinear voltage-current characteristics found in ex- 
periments Q. 

Randomly placed impurities in Ising ferromagnets may 
generate either a random magnetic field or random ex- 
change couplings ||. The free energy landscape for do- 
main walls in these systems is determined by the com- 
petition between the pinning energy and the energy cost 
(per unit length or area) for creating the wall. When the 
system is quenched to a low temperature, the magnetic 
domains grow. As for the flux line, the free energy gain 
due to the motion of a domain wall segment is expected 
to be of the same order as the free energy barrier which 
has to be overcome. 

Since energy barriers play an important role in the dy- 
namics of glassy systems, it is essential to know their 
properties. The scale of these barriers should grow with 
observation size L like a power law . Usually, it is 



assumed that the only energy scale in the system is set 
by the fluctuations in free energy which increase as L e , 
and that therefore ip — 9 However, it is also quite 

possible that the heights of the ridges in the random en- 
ergy landscape scale differently from those of the valleys 
that they separate, with ip > 9. Yet another scenario is 
that transport occurs mainly along a percolating channel 
of exceptionally low energy valleys with ip < 9. 

A first attempt to clarify this situation was taken in 
Ref. p0| , where ip = 9 was established for a FL mov- 
ing in 2 dimensions. Using a combination of analytic 
arguments and numerical simulations, lower and upper 
bounds to the barrier were found. This argument was 
then extended to a FL in 3 dimensions JTTJ , yielding again 
ip = 9. In this article, we present in more detail the ar- 
guments discussed briefly in these earlier papers, includ- 
ing also systems with long-range correlated randomness 
(random-field Ising models) in 2 dimensions. We obtain 
in all cases lower and upper bounds to the barrier that 
scale as L e , except for possible logarithmic factors, lead- 
ing to ip = 9. Furthermore, it is argued that the result 
ip = 9 holds also in higher dimensions, as long as the 
distribution of minimal energies decays exponentially. In 
all cases, the line can move through the system by en- 
countering energy fluctations of only order L s around the 
mean minimal energy. We also show that a line which ini- 
tially has a larger energy can reach this region of minimal 
energies by crossing barriers of order L e (or smaller). 

The outline of the paper is as follows: In section [0], we 
determine the energy barrier for a FL moving in 2 dimen- 



sions. In section [II 



we apply the same algorithm to de- 
termine the energy barrier to the motion of domain walls 
in 2-dimcnsional random-field Ising systems. In section 
|y|, we study energy barriers for a FL in 3 dimensions and 
discuss also the behavior in higher dimensions. In sec- 
tion [v| we take a general look at the energy landscape 
and show that a line can move from any initial configura- 
tion to a minimal configuration by going over no barrier 
higher than L . In section VI, we try to put the defini- 



tion of energy barriers on a more solid foundation, and 
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section VII argues that the results of the paper can be 



generalized to other elastic media with impurities. 



II. ENERGY BARRIERS FOR FLUX LINES IN 2 
DIMENSIONS 



In two dimensions, we represent a DPRM by the fol- 
lowing model: The line is discretized to lie on the bonds 
of a square lattice, directed along its diagonal. Each 
segment of the line can proceed along one of two direc- 
tions, leading to a total of 2* configurations after t steps. 
These configurations are labelled by the set of integers 
{x(t)} for r = 0, 1, • • ■ ,t, giving the transverse coordi- 
nate of the line at each step (clearly constrained such 
that x(t + 1) = x(t) ±1). To each bond on the lattice is 
assigned a (quenched) random energy equally distributed 
between and 1. The energy of each configuration is the 
sum of all random bond energies on the line. Without 
loss of generality, we set x(0) — 0. 

Some exact results are known for this model Qj: The 
fluctuations in the free energy at finite temperature scale 
as i 1 / 3 . The meanderings of the transverse coordinate 
of the line scale as fi, where ( — 2/3 is the roughness 
exponent. The scaling behavior of the pinned FL is gov- 
erned by a zero-temperature fixed point JoJ where energy 
fluctuations scale in the same way. A FL at low tem- 
peratures, and in thermal equilibrium, is likely to spend 
most of the time in configurations of minimal energy. For 
each endpoint (i, x) with x = —t, — t + 2, • • • , t, there is a 
configuration of minimal energy E m i n {x\t) which can be 
obtained numerically in a time of order t 2 . It is known 
that for | a; | < x c oc t 2 / 3 , the function E m i n (x\t) behaves 
as a random walk and is thus asymptotically Gaussian 
distributed flp^|. Since beyond the interval [— x c , x c ] the 
energy of minimal paths is systematically larger, we con- 
sider in this paper only the region [— x c , x c }. Fig.|] shows 
minimal paths of length t = 256 to endpoints between 
x = —96 and x = +96. 

We want to find the energy barrier that has to be over- 
come when the line is moved from an initial minimal en- 
ergy configuration {xi(r)} between (0, 0) and (t, —Xf) to 
a final configuration {xf(r)} between (0, 0) and (t, +£/), 
with Xf = Xf(t) < x c - The only elementary move al- 
lowed is flipping a kink along the line from one side to 
the other (except at the end point). Thus the point (r, x) 
can be shifted to (r, x± 2). Each route from the initial to 
the final configuration is obtained by a sequence of such 
elementary moves. For each sequence, there is an inter- 
mediate configuration of maximum energy, and a barrier 
which is the difference between this maximum and the 
initial energy. In a system at temperature T, the proba- 
bility that the FL chooses a sequence which crosses a bar- 
rier of height Eb is proportional to exp(— Eg/T), multi- 
plied by the number of such sequences. We assume that, 
as is the case for the equilibrium DPRM, the "entropic" 
factor of the number of paths does not modify scaling 



behavior. Thus at sufficiently low temperatures, the FL 
chooses the optimal sequence which has to overcome the 
least energy, and the overall barrier is the minimum of 
barrier energies of all sequences. 

Since the number of elementary moves scales roughly 
as the area between the initial and final lines, the num- 
ber of possible sequences grows as t xt . This exponential 
growth makes it practically impossible to find the barrier 
by examining all possible sequences, hampering a system- 
atic examination of barrier energies. Rather than finding 
the true barrier energy, we proceed by placing upper and 
lower bounds on it. 

The lower bound was given in Ref. |l3], and is ob- 
tained as follows: The endpoint of the path has to visit 
all sites (t,x) with |jc| < Xf, and the energy of any 
path ending at (t,x) is at least as large as E m i n (x\t). 
Therefore the barrier energy cannot be smaller than 
max.[E min (x\t) - E min (-Xf\t)) for x € [—Xf,Xf]. When 
Xf is sufficiently small, the probability distribution of this 
lower bound is identical to that of the maximal deviation 
of a random walk of length x / (lj] . The latter is a Gaus- 
sian distribution with a mean value proportional to y/xj, 
and a variance scaling as Xf. This growth saturates for 
Xf of the order of t 2 / 3 , leading to the scaling behaviors, 



E<?\t,x)) =ff[ sr \x/t i l% 
vax(£i sr) )=t 2 /Vi Sr) (*A 2/3 ), 



and 



(1) 



for the lower bound and its variance. The functions 
f\ ST \y) and f% \V) are proportional to Jy and y for 
small y, respectively, and go to a constant for y = O(I). 
Our simulation results for systems with t = 256, 512, 
1024, 2048, and 4096 confirm this expectation. Fig.| 

shows the scaling functions fi (y) and (y) f° r dif- 
ferent t, and the collapse is quite satisfactory. However, 
the initial growth proportional to y/xj, is not clearly seen 
at these sizes. 

To obtain an upper bound for the barrier, we specify an 
explicit algorithm for moving the line from its initial to fi- 
nal configuration. This is achieved by finding a sequence 
of intermediate steps. It is certainly advantageous to 
keep the intermediate paths as close to minimal configu- 
rations as possible. We therefore proceed in the following 
way: We first find the minimal paths connecting (0, 0) to 
the points (t— 1, x) with Xi(t — I) < x < Xf(t—1), and we 
add a last step to the left (from (t — l,x) to (t, x— 1)). If 
Xi(t— 1) > Xi(t) = —Xf, we then move the point (t, —Xf) 
to (t, —Xf + 2). Now the path has the same endpoint as 
the first intermediate minimal path. We then move the 
path to this first intermediate configuration (the precise 
prescription will be given below), and then we move again 
the endpoint. This procedure is repeated, until the path 
reaches its final configuration. At each step, we obtain a 
local barrier path which separates two neighboring min- 
imal configurations. The overall barrier is of course the 
one with the highest energy. While it may occasionally 
be possible to go to the next intermediate configuration 



2 



in a single elementary move (as defined above), this is 
generally not the case. Intermediate minimal paths with 
the same endpoint may be quite far apart at coordinates 
t < t. The reason is simple: suppose the random poten- 
tial has a large positive fluctuation, a "mountain." The 
minimal energy paths will then circumvent this region by 
going to its right or left. The last path going to the left 
and the first one going to the right have almost the same 
energy. They form a loop which can be quite large and is 
likely to enclose the barrier when both paths separate al- 
ready at small r. Such loops have been conjectured |||7]] 
to play an important role in the low-temperature dynam- 
ics of the DPRM. Since the transverse fluctuations of a 
minimal path of length t grow as t 2 / 3 , we expect the lat- 
eral size of these loops to also be of this order. 

The algorithm for moving a line of length t = 2™ 
from an intermediate configuration {xi(r)} to another 
one {x2(t)}, with x%{t) — X\(t) is as follows: If x% (r) < 
xi(t) + 2 for all r, we can choose a sequence of elementary 
moves such that at most two bonds of the line are not 
on one or the other minimal path, leading to a barrier of 
order 1 between the two. If x 2 (r) > xi(t) + 2 for some 
t, the two paths enclose a loop. We then consider the 
points (t/2 - l,x) withxi(t/2-l) < x < x 2 (t/2-l). For 
each of these points, we find a minimal segment of length 
t/2 — 1, connecting the point (t/2 — l,x) to (0,0) by a 
minimal path, and we take a final step to the left from 
(t/2 — l,x) to (t/2,x — 1). In the same way, we connect 
the points (t/2 + l,x) with x 1 (t/2 + l) < x < x 2 (t/2 + l) 
to x\(t)) via minimal paths and add a first step to 
the right from (t/2,x — 1) to (t/2, x). Two such seg- 
ments form together an almost minimal path of length 
t, constrained to go through intermediate points at t/2 
and t/2 ± 1. We next move the line {xi(t)} step- 
wise through this sequence of almost minimal paths. If 
Xi(t/2 — 1) = Xi(t/2) — 1, we first move the upper seg- 
ment. If Xi(t/2) = Xj(£/2 + l) + l, we then move the lower 
segment. Then we move the middlepoint. We continue 
by repeatedly moving the upper segment, the lower seg- 
ment, and the middle point, until the final configuration 
{x2(t)} is reached. (If the length of the line is different 
from 2 n , we might have to choose the upper segment to 
be larger by 1 than the lower segment, or vice versa.) 

The prescription for moving the segments of length t/2 
is exactly the same as for paths of length t: If the dis- 
tance between two consecutive configurations is larger 
than 2 for some r £ [0,t/2], we consider the points at 
(t/4± l,x) in between the two, and find minimal paths 
of length t/4 — 1 connecting them to the initial and final 
points, and add a step to the middle points. Next we 
attempt to move segments of length t/2 by repeatedly 
moving line portions of length f/4. In some cases, when 
the energy barrier is large, it is necessary to proceed with 
this construction until the cutoff scale of t/2™ -1 = 2 is 
reached. Thus, at each intermediate configuration, the 
line is composed of one segment of length t/2, one of 
length t/4, etc; ending with two smallest pieces of length 
t/2 m (equal to 2 in the worst case). The barrier path is 



the intermediate configuration with highest energy. Fig.|l| 
shows the barrier paths resulting from the above con- 
struction. 

We now estimate the barrier energy resulting from the 
above construction. Each intermediate path is composed 
of segments of minimal paths with constrained endpoints, 
and we would like to find the probability distribution for 
the highest energy. Constraining the endpoints of a min- 
imal path of length r typically increases its energy by 
E^ r \r) oc r 1 / 3 . A subset of these intermediate paths 
(those that cross the largest mountains) have constraints 
imposed on segments of length t, t/2, t/4, and all the 
way down to unity. The number of paths in this subset 
(henceforth referred to as candidate barriers) grows as 
N c (t) oc t a , with 1 < a < 1 + 2/3. The lower limit comes 
from noting that for each loop of size 2 m there exist at 
least two loops of size 2" i ~ 1 , one in the upper and one 
in the lower half of the parent loop, thus N c > t. The 
upper limit comes from the total number of intermedi- 
ate configurations that grows as txf. The energy of each 
candidate barrier path is obtained in a manner similar to 
that of the lower bound: Instead of finding the maximum 
of a random walk of length xt oc t 2 / 3 , we now have to 
examine the sum of the maxima for a sequence of shorter 
and shorter random walks added together. The mean 
value of this sum is related to the convergent series, 

(4 sr) (t)) = 

= (E ( _ sr) (t) + 2E ( _ sr) (t/2) + 2E { _ sr) (t/A) + ■ ■ A + Aln(t) 
= (^i sr) (t)) (l + 2(2- 1 / 3 + 2- 2 / 3 + ---)) +Ahi(t)+B 
~ (£i sr) (t)) (-l + 2(l-2- 1 / 3 )- 1 ) + A\u(t)+B 
= 8.m...(E {sr) (t)) +A\n(t)+B. (2) 

The correction term Aln(t), is explained as follows: Each 
segment of length 2™ is composed of a minimal path of 
length 2 rn — 1 and one step which has a random energy 
(the final or initial step, depending on whether the seg- 
ment lies in the upper or lower half of a loop). So the 

energy of the segment is equal to the energy E^ r '(2 m ) of 
a minimal path of length 2 m , plus a constant of order 1. 
Since a candidate barrier has n — ln(t) / ln(2) segments, 
these constants add up to A ln(t) , with A of the order of 
unity. The constant B in Eq.(g) accounts for the break- 
down of the scaling form of the energy increase for small 
loops. The mean angle of the smallest loops (of size 2) 
approaches the 45° limit; their mean energy growing as 
0.5t m . For the larger loops, the angle tU 3 /t m is small and 
the energy is 0.23t m . A finite value of m acts as a cutoff 
separating the two limits. The energy difference per unit 
length between small and large paths then leads to the 
additive constant B (of the order of unity) in Eq. (|J) . 

The barrier energy is the maximum of the N c (t) ener- 
gies of all candidate barriers. To find its characteristics, 
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we need the whole probability distribution for the en- 
ergy Ec Sr \t). Since Ec is the sum of energies coming 
from its minimal segments, the simplest assumption is 
to regard the segment energies as independent, approx- 
imately Gaussian, random variables. We then conclude 
that i?c Sr ' (t) is also Gaussian distributed with a variance, 

var (^)(t)) 

= var (^i sr) (<)) + 2 (var (V sr) (t/2)) + • • ■) 

~ 4.40... var (^i sr) (t)) ott 2 / 3 . (3) 

Since the different segments are in fact constructed 
through a specific recursive procedure, their indepen- 
dence cannot be justified. In the worst case that they 
are completely dependent, the right-hand side of Eq.(||) 
has to be multiplied by n = log 2 (i). Since our numerical 
results show no evidence for such a logarithmic factor, 
we shall not consider it any further. 

It can be checked easily that (for large TV), the max- 
imum of N independent Gaussian variables of mean a 
and variance a 2 , is a Gaussian of mean a + a\/2\aN and 
variance a 2 /(2lnN) |l5[ . Since the N c (t) candidate bar- 
riers have large segments in common, their energies are 
not independent. We can approximately take this into 
account by assuming a subset of them as independent, 
leading to N cx t a for some a' < a. We thus obtain the 
following estimates for the mean upper bound in barrier 
energy, 

(£<r>(ar,t)) = (Ejr\x,t)) + yj2lnNvaxEi sr \x,t) 
^(l + /3^V^l)t^g[ sr \x/^% (4) 
and its variance, 

var(.El sr) 0M)) i r -^- 

^atrHx/&*). (5) 

The functions g[ sr \y) and g^ r \y) are proportional to 
y/y and y respectively, for small y, constant at large y, 
and in general different from those of the lower bound. 

Our numerical simulations indeed confirm the above 
scaling forms. The scaling functions g[ sr \y) and g^ r \y) 
are plotted in Fig.|| for different values of t, after aver- 
aging over 2000 realizations of randomness. The y/ln(t) 
factors are essential, as a comparable collapse is not ob- 
tained without them. In fact, the best fit to < E^ r \t) > 
is obtained by including the correction to scaling term 
cx< E^ r '(t) >, and with = 1. The numerics 

therefore support the neglect of correlations, and the as- 
sumption of a Gaussian distributed Ec Sr \t). As in the 
lower bound, the initial scaling proportional to ^Jx] is 



not clearly seen for the sizes studied. Since the leading 
power for the scaling of the lower and upper bounds is 
identical, we conclude that the barrier energies also grow 
as t x l z . (It remains to be seen if the logarithmic factors 
are truly present, or merely an artifact of our algorithm.) 

III. ENERGY BARRIERS FOR DOMAIN WALLS 
IN 2-DIMENSIONAL RANDOM FIELD ISING 
SYSTEMS 

In the previous section, we considered random bond 
energies which were uncorrelated. The analytic argument 
for the upper bound relied on the random-walk behavior 
of E m i n (x\t) in this situation. Thus, the proof for ip = 8 
can not directly be extended to other situations, where 
the distribution of lower bound energies is not known. 
An important example is the case of domain walls in 2- 
dimensional random-field Ising magnets. The energy for 
creating a domain wall is equal to the cost of flipping all 
spins on one side of the interface; in turn proportional to 
the sum of all random fields on the flipped spins. There 
are consequently long-range correlations in the domain 
wall energy in the direction perpendicular to the orien- 
tation of the wall . 

We describe the configurations of the domain wall by 
essentially the same model as the FL, but assigning to 
each bond a random energy with long-range correla- 
tions in the x-direction. These correlations are gener- 
ated by first selecting for each time t, random numbers 
{r\ N \r\ N+1 \ ■ ■ ■ ,r[ N 3 \r[ N ^} equally distributed 
between —1 and 1, where N is (at least) as large as the 
largest time occuring in the simulations. To each bond 
connecting (t, x) to (t+1, x±l) we then assign the energy 

/s-l/2±l/2 itf-l \ 

tM e rt ~ £ r * ■ 

V \ i=— JV x+l/2±l/2 / 

Fig.|| shows minimal paths of length t — 128. Due 
to the correlations, neighboring bonds have almost the 
same energy, and therefore minimal paths tend to have 
large parallel portions. Fig.|| shows the minimal energy 
as function of the endpoint position for a given realiza- 
tion of randomness, and for t — 1024. This curve is 
much smoother and has longer correlations than the cor- 
responding curve in the case of short-range correlated 
randomness, where the minimal energy performs a ran- 
dom walk. 

The fluctuations in free energy of a line are known to 
scale as t, and the roughness exponent is £ = 1 p7| . 
We determined numerically the distribution function for 
the minimal energy shown in Fig.|5|. It is very close to 
a Gaussian, with no apparent power-law tails. We will 
show that, due to this property of the minimal energy 
distribution, the lower and upper bounds to the barrier 
scale in the same way as the fluctuations in minimal en- 
ergy. 
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As in the previous section, we move the line from 
an initial minimal energy configuration {xi(r)} between 
(0,0) and (t, — Xf), to a final configuration {x/(r)} be- 
tween (0,0) and (t,Xf). Since the endpoint of the path 
has to visit all sites (t, x) with \x\ < Xf, and since the 
energy of any path ending at (t, x) is at least as large 
as E m i n (x\t), the barrier energy cannot be smaller than 
m.ax[E m i n (x\t) — E rn i n (—Xf\t)]. Since the distribution of 
minimal energies decays exponentially and has no power- 
law tails, we can expect that the lower bound scales in 
the same way as the fluctuation of the minimal energy, 
leading to 

(E { l r \t,x)\ =tf[ lr \x/t), and 

™(E^) = t 2 ft\x/t), (6) 

for the lower bound and its variance. Our simulation 
results for systems with t = 256, 512, 1024 and 2048 
confirm this expectation. Fig.^J shows the scaling func- 
tions /{ (y) and f^iu) for different t, and the collapse 
is quite satisfactory. However, the initial growth propor- 
tional to Xf, is not clearly seen at these sizes. Fig]?] shows 
the distribution of lower bound energies. It is very close 
to a (half)-Gaussian of width proportional to t. 

An upper bound can be obtained by exactly the same 
algorithm as before. The analytic argument made in the 
previous section, however, cannot be directly repeated, 
since the function E m i n {x\t) is no longer a random walk 
in x, and since we do not have analytic results for the 
lower bound. We can, however, combine analytic argu- 
ments with the numerical results for the lower bound to 
predict the scaling behavior of the upper bound. Since 
the line is always composed of minimal segments, the 
energy of a candidate barrier which has segments of all 
lengths down to the cutoff is given by 

(2e j w) 

= (E { l r) (t) + 2 (E { l r \t/2) + E { l r) (t/A) + •••)) 
~3^i' r) (t)) +A'ln(t)+B'. (7) 

The origin of the terms A' \n(t)+B' has been explained in 
the previous section (see paragraph after Eq.(^|)). Since 
the energy distribution of the lower bound is approx- 
imately (half-)Gaussian, the energy distribution of the 
candidate barriers decays also like a Gaussian. The up- 
per bound to the barrier energy is the maximum of the 
energies of all candidate barriers. In our simulations, we 
find no evidence for logarithmic factors, indicating that 
the number of candidate barriers increases either very 
slowly, or not at all, with t. From Fig.|| we can see that 
there is essentially one large loop over a distance of the 
order of the length of the path, leading to only few inde- 
pendent candidate barriers. As in the previous section, 
we find that the maximum of these candidate barriers, 
which is the upper bound to the barrier energy, scales in 
the same way as the lower bound, i.e. 



(^ r) (x,t)) = (EW(x,tj) + ^21n^(M V arSf \x,t) 
^tg{ lr) (x/t), (8) 

and its variance scales as 

var(j5$ r) (a:,t)) = wfe fr )(x,i)) /21nA^ r > 

*1?g« r \x/t). (9) 

Fig.^ shows the scaling functions f% and g^ ■ 

To summarize the results so far, we have established 
the relation ip — 9 for lines in 2-dimcnsional systems 
with short- and long-range correlated randomness. Since 
the considerations for both systems rely strongly on the 
dimensionality, it is of importance to look also at a 3- 
dimensional system, which is physically more relevant. 



IV. ENERGY BARRIERS FOR A FLUX LINE IN 
3 DIMENSIONS 

In a two-dimensional system, the endpoint of the FL 
has to move through all points (x,t) with Xi < x < Xf. 
This property was essential for the derivation of the lower 
bound in the previous sections. A FL which moves in 
three dimensions can avoid regions in space which are 
energetically unfavorable for some of its segments, and 
one may therefore speculate that ip < 6. In this section, 
we first determine numerically a lower bound for the bar- 
rier energy which scales in the same way as the energy 
fluctuations, thus ruling out ip < 9. Further numerical 
results predict that an upper bound scales in the same 
way, thus leading to ip = 9. 

The line now lies on the bonds of a cubic lattice, start- 
ing at the origin and directed along its (1,1,1) diagonal. 
Each segment of the line can proceed in the positive di- 
rection along one of the three axes, leading to a total 
of 3* configurations after t steps, with endpoints lying in 
the plane which is spanned by the points (t, 0, 0), (0, t, 0), 
and (0,0,i). A given configuration of the FL is labelled 
by vectors {x(t)} for r = 0, 1, • • • , t, giving the transverse 
coordinates of the FL at each step. The points {x(t)} lie 
on the vertices of a triangular lattice. For a given value 
of t, they lie on one of three alternating sublattices. 

The minimal energy E m i n (x\t) can be obtained numer- 
ically in a time of order t 3 . The fluctuations in minimal 
energy are known to scale as t e with 9 ~ 0.24, and the 
roughness exponent for minimal paths is C — 0.62 fll"8| , |l9| ] . 
The endpoints of the minimal paths with the lowest en- 
ergy lie within a distance proportional to from the 
origin. Figure ^ shows the minimal energies of paths of 
length t = 288 to endpoints x with \x\ < 0(t c ). The 
highest energy in this region is represented in white, the 
smallest energy in black. The minimal energies are corre- 
lated over a distance of the order of tf- . The distribution 



5 



of minimal energies is close to a Gaussian and is shown 
in FigJ|. Similar to a 2-dimensional system (see also 
Fig.||), this distribution seems to have a third cumulant 
since it is not completely symmetric. 

A lower bound to the barrier energy is obtained as 
follows: While the line moves from its initial to final 
configuration, the transverse coordinates of its endpoint 
move between nearest-neighbor positions on one of the 
above mentioned triangular sublattices. When the end- 
point is at a position x, the energy of the line is at least 
as large as the minimal energy E m i n (x\t) . The maximum 
of all these minimal energies along the trajectory of the 
endpoint, minus the energy of the initial configuration, 
certainly bounds the barrier energy from below. Since 
we do not know the actual trajectory of the endpoint, we 
have to look for the trajectory with the smallest maximal 
energy. Only in this way can we be sure that we have 
indeed found a lower bound. This situation is fundamen- 
tally different from a 2-dimensional system, where there 
is only one possible trajectory for the endpoint. 

Provided that the minimal energies E m i n (x\t) are 
known, this lower bound is determined in polynomial 
time by using a transfer-matrix method: We start by as- 
signing to the initial point Xi a "barrier energy" B(xi) = 
0, and to all other sites x on the same sublattice a bar- 
rier energy B(x) = t, which is certainly larger than the 
lower bound resulting from the algorithm after many it- 
erations. At each step, the energy B of all sites x, except 
for the initial site, is updated according to the following 
rule: Look for the minimum of the energies B(x ± e*;) 
of the 6 neighbors. If this is smaller than B(x), replace 
B(x) by this minimum, or by E min (x\t) - E min (x l \t), 
whichever is larger. After a sufficiently large number of 
iterations, which is of the order of the size of the area of 
interest (scaling as i 2 ^), all possible trajectories to end- 
points within this area have been probed, and the barrier 
energies B(x) do not change any more. The energy B{xf) 
is then identified as the lower bound. Figure shows the 
lower bound to the energy barrier for a line with the end- 
point moving from the origin to sites within a distance of 
the order of ft, for different values of t and averaged over 
500 realizations of randomness. The distance \xf — Xi\ 
has been scaled by ft, and the energy by t e . After this 
rescaling, all the curves collapse, leading to the following 
scaling behavior for the lower bound, 

(£_(*, \x f - = t 9 f-(\x f - Xi\lt<). (10) 

The function f(y) should be proportional to y e ^ for 
small y. Again, for the simulated system sizes, this 
asymptotic scaling is not clearly seen. For y > 1, the 
scaling form in Eq.([l0|) breaks down since the minimal 
energy is then a function of the angle (|ir|/i). We con- 
clude that the lower bound to the barrier scales in the 
same way as the fluctuations in minimal energy, and con- 
sequently the energy barrier increases at least as t , lead- 
ing to tp > 6. The distribution P(E_) of the lower bound 
energy for a fixed distance \x\ cx ft is shown in Fig.jn]. It 



appears to be half-Gaussian with width proportional to 
ft. 

The result %p > 6 is not surprising if we note that an 
even simpler lower bound is given by max(E m i n (xf\t) — 
E m in(xi\t),0), which evidently scales as t e since the dis- 
tribution function of minimal energies decays exponen- 
tially fast, i.e. has no power-law tails (see Fig|j). To 
make sure that the scaling of the lower bound found 
above is not dominated by the neighborhood of final con- 
figurations with particularly high energies, we repeated 
the above simulations by allowing only endpoints with 
minimal energies smaller than the initial energy. This 
corresponds to a situation where the endpoint of the line 
only moves to positions which are energetically more fa- 
vorable. The results are shown in FigJlO] and again col- 
lapsed by the scaling form 

(E-(t, \x f - Xi\)) = t 9 f-(\x f - Xi\/^). (11) 

As in the previous case, the asymptotic scaling f-(y) oc 
y s /C f or small y cannot be clearly seen. The energy dis- 
tribution of the lower bound is again a half-Gaussian of 
width proportional to t^> and looks similar to Fig.[TT|. 

The same scaling behavior is also found when instead 
of the optimal trajectory for the endpoint, the shortest 
trajectory (a straight line) is chosen. In this case, the 
mean of the barrier energy Eq has the scaling form 

(Eo(t, \x f - Si\)) = t e f (\x f - Xi\/t<) (12) 

(see Fig.|l0|) , again with a half-Gaussian distribution of 
width proportional to ft . This, of course, does not rep- 
resent a lower bound to the true barrier, but it will be 
important for the determination of an upper bound be- 
low, and is therefore included here. 

The result E_ cx t 6 (Eq.([ll])) can be explained from 
the exponential tails of the distribution of minimal en- 
ergies: If we asume that the endpoint of the line moves 
only in valleys of particularly low energy, we can succes- 
sively remove all sites with the largest minimal energy 
from the set of possible endpoints, until the connectiv- 
ity over the distance ft breaks down. The remaining 
endpoints form percolation clusters, and their density is 
given by the corresponding percolation threshold. (This 
is analogous to random resistor networks describing the 
hopping resistivity for strongly localized electrons. The 
resistance of the whole sample is governed by the critical 
resistor that makes the network percolate plf .) Since 
the occupied sites are correlated over the distances con- 
sidered, the value for the threshold is different from the 
site percolation limit of 0.5 in an infinite triangular lat- 
tice with no correlation between occupied sites. But for 
the present purpose, it is sufficient to know that this 
threshold is finite, and that therefore a finite percentage 
of all sites are below threshold. Since the distribution 
of minimal energies decays rapidly, its tail cannot con- 
tain a finite percentage of all sites. We conclude that the 
threshold is within a distance of t from the peak, and 
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therefore that the energy fluctuation on the percolation 
cluster, and consequently the lower bound for the barrier, 
are proportional to t e . 

We now proceed to construct an upper bound to the 
energy barrier. To this purpose, we specify a sequence of 
elementary moves which take the line from its initial to 
final configuration. The only elementary move allowed 
is flipping a kink along the line. Thus the point (r, x) 
can be shifted to (r, x ± e^), where ±e*j are the six vec- 
tors which connect a vertex in the triangular lattice to 
its nearest neighbors within the same sublattice. The al- 
gorithm is similar to the one in 2 dimensions: First, we 
choose a sequence of endpoints connecting the initial to 
the final endpoint which is as short as possible. Then, we 
draw all the minimal paths leading to these endpoints, 
and attempt to move the line through them sequentially. 
If two consecutive minimal paths have nowhere a distance 
larger than 1 (measured in units of | | ) , we can choose 
a sequence of elementary moves such that at most two 
bonds of the line are not on one or the other minimal 
path, leading to a barrier of order 1 between the two. If 
the distance is larger than 1, we proceed essentially in the 
same way as in 2 dimensions, i.e. we consider the midway 
points (t/2,Xi) which connect both lines via the shortest 
possible trajectory Xi(t/2) (if there are several possibili- 
ties, we choose one at random). For each of these points, 
we find two minimal segments of length t/2 connecting on 
one side to (0,a;) and on the other to either (t,X\(t)) or 
(t,X2(t)). Then wc move the line by repeatedly moving 
segments of length t/2, etc. 

The energy of a candidate barrier is then given by 



(E ( - 3d ' ) (t,t c )) 

^(£ (M C ))(l + 2((l/2) e + (l/4) e 
= (£o(M C )) (-l + 2/(l-(l/2) e )) 
- 12.0... (£ (M C )>. 



(13) 



In principle, one should add correction terms similar to 
those in EqsJ| and (t| However, these corrections are sub- 
leading with respect to t e and will be neglected. 



As mentioned in section HI, we cannot rule out that 



the energies of minimal segments are independent from 
each other. In the worst case, where they are completely 
dependent, Eq.(|l3|) has to be multiplied by \/log 2 t. This 
may result in an additional factor proportional to Vint 
in the upper bound, but does not otherwise affect any 
of our conclusions. The number of independent candi- 
date barriers increases with some power in t. Since their 
energy distribution decays like a Gaussian, we can take 
their maximum in the same way as before, and we finally 
obtain the following estimate for the upper bound, 



E^\t,\x f 



Xi\) 



= (E c (\x f -Xi\,t)) + V2hTNva,rE c (\x f 

(Vhat)^/+(I^/-^IA C )- 



(14) 



We have thus shown that the energy barrier encoun- 
tered by a FL moving in a 2d or 3d random medium 
has an upper and a lower bound which both increase as 
t , except for logarithmic factors. It thus follows that 
the barrier itself scales as t e , confirming the hypothesis 
tp = 9. Since the arguments are mainly based on the 
exponential tails of the minimal energy distributions, it 
is expected that the result ip — 8 holds also in higher 
dimensions. The only requirement is that the tails in the 
distributions of minimal energies still decay sufficiently 
rapidly. 



V. BARRIERS TO FAR FROM MINIMAL 
CONFIGURATIONS 

In the previous sections, we discussed energy barriers 
which have to be overcome by a line moving between min- 
imal energy configurations. We showed that such lines 
can stay in an energy interval (E m i n ) ± const t e . How- 
ever, a line may initially have an energy which is much 
larger. The initial configuration of a FL penetrating the 
system may be straight and parallel to the external mag- 
netic field. If a system is cooled down from high tem- 
peratures, configurations of the FL are random walks of 
roughness exponent ( = 1/2. An initial configuration 
with roughness exponent £ = 1 is found for FLs driven 
close to a depinning transition [^2| . If the temperature 
is low (as we always assume in this paper), the line then 
relaxes to some metastable state. The FL will ultimately 
reach a configuration of minimal energy, only if it is not 
separated from it by abnormally high barriers. We there- 
fore show in this section that the line can reach the min- 
imal energy region by going only over barriers which are 
not larger than order of t e . We specify an algorithm for 
moving a line of length t = 2 n from any initial configura- 
tion to one of minimal energy. The algorithm is similar 
to that presented in the previous sections, and leads to 
barriers of the order of t . 

First, we assume that its initial roughness is not larger 
than that of minimal energy paths. Let {x u (t)} for 
t = 0, • • • , t be the initial configuration of the line, and 
{xq(t)} a minimal energy configuration with x n (0) = 
xq(0) and x n (t) = xo(t). We then define a sequence of 
paths {a; m (r)}, m = 1, • • • , n— 1, which are constrained to 
go through the points x n (kt/2 m ) for k = 0, 1, • ■ • , 2 m and 
are composed of 2 m minimal segments of length t/2 m . 
The energy of such a segment is smaller than the energy 
of any other piece of a path with larger m which has the 
same endpoints as the segment. We now move the line 
successively through this sequence of configurations, go- 
ing from the largest to the smallest value of m. The con- 
figurations {x m+ i(r)} and {x m (T)} intersect each other 
at the points r = kt/2 m , with k = 0, ■ • • , 2 m . We there- 
fore can move the line from the configuration {a; m +i(r)} 
to {x,„(t)} by successively moving segments of length 
t/2 m . In many cases, the segments have to overcome a 
loop, and then we apply the algorithm defined previously. 
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In contrast to the previous sections, these loops do not 
separate two minimal configurations, but one minimal 
segment, and another constrained at its midpoint, a con- 
stellation which occured also in the previous sections as 
an intermediate situation. Since we restricted the rough- 
ness of the initial configuration to less than that of min- 
imal paths, the size of the loops does not exceed t*. The 
number of independent candidate barriers within a loop 
is therefore smaller than, or equal to, N c oc (2 m ) a from 
previous arguments, where the exponent a' depends on 
the model. The energy of each candidate path is smaller 
than, or equal to, the energy E c (2 m ), which was also 
obtained in the previous sections. The total number of 
loops is less than or equal to 1 + 2 + • • • + 2 ,l ~ 1 < 2™ = t, 
and the energy of each candidate barrier is certainly over- 
estimated if we assume that all loops are of size t. We 
therefore find an upper bound to the barrier which is 
the maximum of tt a — t 1+a candidates chosen from 
a distribution P(E c (t)) with (E c (t)) oc t e , and with a 
Gaussian tail. As we saw in section III, such a maxi- 
mum scales as t 6 



'Int. We therefore have shown that 
a line can move from any configuration with roughness 
exponent less than £ to a minimal energy configuration 
by crossing barriers which are not larger than order of 
t e , provided that the barriers between minimal configu- 
rations scale also as t 6 . 

A similar result can be obtained for any initial con- 
figuration of the line. To demonstrate this, let us look 
at the configuration x t (r) = — r, which is as far as pos- 
sible from a minimal configuration. We then define a 
sequence of paths {x m (r)}, for to = t — 2, t — 4, • • • , 0, 
with x m (r) = —t for r < m, and connecting the points 
(to, —to) and (t, —to) by a minimal path. To go from one 
configuration to the next one, the line has to overcome 
a loop of size no bigger than (i — m) e < t e . There are 
consequently proportional to t 1+a candidates for barrier 
paths of length between 2 and t. We certainly find an 
upper bound to the barrier by assuming that all these 
candidates have the length t, and that their energies are 
taken from a (half-) Gaussian distribution of width pro- 
portional to t s . The upper bound consequently scales as 
t e VhTt. 



VI. MULTIPLE BARRIERS 

So far, we tacitly assumed that the activation barrier is 
given by the difference of the highest energy encountered 
by the line and its initial energy, just as in thermally acti- 
vated chemical reactions. This assumption, however, has 
no solid foundation, since the line does not simply move 
over an isolated maximum, but through a random energy 
landscape. In addition, it is not at all clear how results 
obtained for a point-like particle in a 1-dimensional en- 
ergy landscape can be generalized to lines moving in 2- 
or 3-dimensional systems. To shed some light at least 
on the first of these points, we study in this section a 



particle in a 1-dimensional energy landscape at low tem- 
peratures. We tilt this landscape by a small angle to take 
into account the effect of an external driving force. Using 
the Fokker-Planck equation, we calculate the stationary 
particle current through this tilted energy landscape. We 
find that it is not the difference between the maximal and 
initial energies, but the difference between the maximal 
and the minimal energies, which determines the activa- 
tion barrier. 

The Fokker-Planck equation for the probability den- 
sity P(x,t) of an overdamped particle in one dimension 
is 



dP _ r _9_ 

dt dx 



kT 



d dV(x) 



dx 



dx 



P(x,t). 



(15) 



Here, T is the inverse of the product of the particle mass 
and the friction coefficient. The potential V(x) is the 
sum of the random potential Vb (x) and a driving term 
—Fx, where F is the constant driving force. Depending 
on the boundary conditions, this equation has different 
stationary solutions dtP = 0. If the boundary is an in- 
finitely high wall at both ends of the system, we obtain 
the equilibrium solution P(x) oc exp(— V(x)/kT), and 
zero current 



i = r 



dx 



dV(x) 
dx 



P{x,t) =0. 



We instead look for a solution where particles enter the 
system at one end and leave it at the other end. This solu- 
tion is most readily found by assuming periodic boundary 
conditions, P(L) = P(0) and V B (L) = V B (0). This situ- 
ation corresponds to a periodic energy landscape which 
has been tilted, and where each section of length L con- 
tains the same number of particles. Clearly, this leads to 
a stationary flow through the system, with particles en- 
tering a section at one end and leaving it at the other. In 
the limit of small driving force F, the stationary current 
is found by considering only terms up to linear order in F 
(order zero gives the equilibrium solution of the untilted 
system), and is given by [B3| 



j = FFL/ 



L e~ v -/ kT dx f L e +VB ' kT dx 



(16) 



In the limit T — > 0, the integrals are dominated by the 
neighborhoods of the maximum and the minimum of the 
potential, leading to j oc exp[(Vg lax - V^ nm )/kT}. This 
means that the particle mobility is determined by the 
difference between the energy maximum and minimum. 
This result is plausible since the particles explore all of 
the energy landscape and therefore also go down to the 
valleys and have to come up all the way again Q . 

Generalizing the above arguments to a line in 2 dimen- 
sions is difficult, and we did not succeed in solving the 
corresponding Fokker-Planck equation analytically. Evi- 
dently, the line can avoid configurations with particularly 
high energy, which seems to justify the assumption that 
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the barrier is the lowest possible which separates the ini- 
tial and the final configurations. In the light of Eq.(16), 
however, we may need to define the barrier energy as 
the difference between the maximum and the minimum, 
instead of the difference between the maximum and the 
initial energy. If so, we should add to the barrier the dif- 
ference between the initial energy and that of the abso- 
lute minimum along the trajectory of the line. We know, 
however, that the distribution of minimal energies has 
only exponental tails, and that therefore both types of 
barriers scale in the same way. Consequently, our results 
do not depend on the precise definition of the barrier. To 
confirm this, we plot in Fig.JTJjthe scaling functions g[ sr ^ 

-(lr) 

and g\ for the barrier to 2-dimcnsional lines with cither 
short-range or long-range correlated randomness, deter- 
mining the difference between the maximal and minimal 
energy of all intermediate configurations of the line. It 
can clearly be seen that the barrier energy still satisfies 
Eqs.(g) and (§). 



VII. CONCLUSIONS 

In this paper we considered various properties of the 
energy landscape of one of the simplest realizations of 
glassy systems. We showed that, under fairly general 
conditions, the energy barriers encountered by a line de- 
scending into the region of minimal energies, or moving 
within this region, scale in the same way as the fluctua- 
tions in minimal energy. This means, in particular, that 
there exist no metastable configurations which cannot be 
left by going over energy barriers smaller than, or equal 
to, the fluctuations in minimal energy. 

Similar arguments are applicable to interfaces in ran- 
dom media, like domain walls in 3-dimensional random 
bond and random field Ising models. When the interface 
moves from an initial to a final configuration, with part of 
its boundary fixed, a given boundary point moves along 
a line. For each position of this boundary point, there 
exists a configuration of minimal energy. The maximum 
of these minimal energies certainly is a lower bound to 
the barrier. If the distribution of minimal energies has no 
power-law tails, this lower bound scales in the same way 
as the minimal energy fluctuations. An upper bound can 
be constructed using a similar algorithm as for the line: 
Each time the interface (or a segment of it) has to over- 
come a loop, we bisect it and repeatedly move the upper 
and lower segment through a sequence of minimal con- 
figurations. In this way, the interface is always composed 
of minimal segments, and it should scale in the same way 
as the lower bound (except for logarithmic factors, and 
provided that the lower bound energy distribution has no 
power-law tails). 

Given these results for lines and interfaces, it is likely 
that they generally hold for elastic media with impuri- 
ties, e.g. for a bunch of flux lines. The latter situation 
is certainly of much more physical releance than a single 



FL. Our results for a single FL, may thus have provided 
a glimpse into the complexity of the energy landscape of 
more complicated glassy systems. 

Based on results for particles in 1-dimensional energy 
landscapes, we also argue that the energy barrier should 
not be defined with respect to the initial energy, but to 
the minimal energy along the trajectory of the line. It 
stills remains a challenge to generalize this argument to 
lines in 2- or higher dimensional energy landscapes, and 
to find a more precise expression for the response of the 
line to a driving force, starting from the Fokker-Planck 
equation. 
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FIG. 1. Minimal paths of length t = 256 for a FL in 2 di- 
mensions to endpoints between x = —96 and x = +96 (solid), 
and the barrier paths between them (dotted). 
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FIG. 2. Scaling functions f[ ar \y), ft\y), O^iv), and 
g^ r \y) (see Eqs.(l), (4), and (5)) for the mean and variance 
of the lower and upper bounds; averaged over 2000 realiza- 
tions of randomness, for t = 256 (solid), 512 (dotted), 1024 
(dashed), 2048 (long dashed), and 4096 (dot-dashed). 




FIG. 3. Minimal paths of length t — 128 to endpoints 
between x — — t and x = +t (solid), and the barrier paths 
between them (dotted), in the random-field Ising model. 
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FIG. 4. The minimal energy E m i n as function of the end- 
point position x for t = 1024 and (a) long-range correlated 
randomness, (b) short-range correlated randomness. 
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FIG. 5. The distribution P(E m i n ) of minimal energies for 
t = 1024 and long-range correlated randomness. The solid 
line is a Gaussian distribution; deviations from it indicate a 
third cumulant. 
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FIG. 6. Scaling functions f[ lr \y), ft\y), g[ lr) (y), and 
92 T \v) ( see Eqs.(6), (8), and (9)) for the mean and variance 
of the lower and upper bounds; averaged over 1000 realiza- 
tions of randomness, for t = 256 (solid), £512 (dotted), 1024 
(dashed), and 2048 (long dashed). The straight line has the 
slope 1. 



FIG. 7. Distribution of minimal energy for an Ising do- 
main wall in 2 dimensions with random fields. The symbols 
are the same as in the previous figure. 




FIG. 8. Minimal energies of paths of length t = 288 in 3 
dimensions to endpoints x with \x\ < O^). The higher en- 
ergies are indicated by white shading, and the lower energies 
by darker shades. 
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FIG. 9. Probability distribution P(E min ) of minimal ener- 
gies E m i n (0\144) in 3 dimensions, averaged over 50,000 real- 
izations of randomness. The solid line is a Gaussian distribu- 
tion. 




FIG. 10. Scaling functions f-{y), f-(y), and fo(y) de- 
fined in Eqs. (10)-(12) for t = 72 (solid), 144 (dotted), 288 
(dashed), 576 (long dashed), and 1152 (dot-dashed), aver- 
aged over 500 realizations of randomness. The straight line 
has slope 0/£ = 0.39. 
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FIG. 11. Probability distribution of The parameters 
and symbols are the same as in Fig. 10. 




10' 3 10" 2 10'' 10° 10' 



FIG. 12. Scaling functions g\ (y) and g\ (y), similar to 
g[ sr \y) in eq.(4) and g[ lr \y) in eq.(8), but with the barrier 
definition of section VI. 
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